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We obtain, by extensive direct numerical simulations, time-dependent and equal-time structure 
functions for the vorticity, in both quasi-Lagrangian and Eulerian frames, for the direct-cascade 
regime in two-dimensional fluid turbulence with air-drag-induced friction. We show that differ- 
ent ways of extracting time scales from these time-dependent structure functions lead to different 
dynamic- multiscaling exponents, which are related to equal-time multiscaling exponents by different 
classes of bridge relations; for a representative value of the friction we verify that, given our error 
bars, these bridge relations hold. 
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The scaling properties of both equal-time and time- 
dependent correlation functions close to a critical point, 
say in a spin system, have been understood well for 
nearly four decades [241125]. By contrast, the develop- 
ment of a similar understanding of the multiscaling prop- 
erties of equal-time and time-dependent structure func- 
tions in the inertial range in fluid turbulence still remains 
a major challenge for it requires interdisciplinary studies 
that must use ideas both from nonequilibrium statistical 
mechanics and turbulence [26-35 . We develop here a 
complete characterization of the rich multiscaling prop- 
erties of time-dependent vorticity structure functions for 
the direct-cascade regime of two-dimensional (2D) tur- 
bulence in fluid films with friction, which we study via 
a direct numerical simulation (DNS). Such a characteri- 
zation has not been possible hitherto because it requires 
very long temporal averaging to obtain good statistics 
for quasi-Lagrangian structure functions [36 , which are 
considerably more complicated than their conventional, 
Eulerian counterparts as we show below. Our DNS study 
yields a variety of interesting results that we summa- 
rize informally before providing technical details and pre- 
cise definitions: (a) We calculate equal-time and time- 
dependent vorticity structure functions in Eulerian and 
quasi-Lagrangian frames [36]. (b) We then show how to 
extract an infinite number of different time scales from 
such time-dependent structure functions, (c) Next we 
present generalizations of the dynamic-scaling Ansatz, 
first used in the context of critical phenomena [25] to 
relate a diverging relaxation time r to a diverging corre- 
lation length ^ via r ~ , where z is the dynamic-scaling 
exponent. These generalizations yield, in turn, an infin- 
ity of dynamic- multiscaling exponents [28] [29] [3Tti35] . (d) 
A suitable extension of the multifractal formalism [27] . 
which provides a rationalization of the multiscaling of 
equal-time structure functions in turbulence, yields lin- 



ear bridge relations between dynamic-multiscaling expo- 
nents and their equal-time counterparts [28l |29l [3TH35] ; 
our study provides numerical evidence in support of such 
bridge relations. 

The statistical properties of fully developed, homoge- 
neous, isotropic turbulence are characterized, inter alia, 
by the equal-time, order-p, longitudinal-velocity struc- 
ture function Sp{r) = {[Su\\{r,t)]P) , where Su\\{r,t) = 
[u(x + r,t) — u(x, t) • r/r], u(x, t) is the Eulerian ve- 
locity at point X and time t, and r =| r |. In the 
inertial range t^^^ <C r <C Sp{r) ^ r^p , where Cp? 
r]d^ and L, are, respectively, the equal-time exponent, 
the dissipation scale, and the forcing scale. The pio- 
neering work [26] of Kolmogorov (K41) predicts simple 
scaling with C^^^ = p/3 for three-dimensional (3D) ho- 
mogeneous, isotropic fluid turbulence. However, experi- 
ments and numerical simulations show marked deviations 
from K41 scaling, especially for p > 4, with (p a non- 
linear, convex function of p; thus, we have multiscaling 
of equal-time velocity structure functions. To examine 
dynamic multiscaling, we must obtain the order-p, time- 
dependent structure functions J'p(r, t), which we define 
precisely below, extract from these the time scales rp(r), 
and thence the dynamic-multiscaling exponents Zp via 
dynamic-multiscaling Ansatze like rp{r) ^ r^p. This task 
is considerably more complicated than its analog for the 
determination of the equal-time multiscaling exponents 
(p [28 - 35 for the following two reasons: (I) In the conven- 
tional Eulerian description, the sweeping effect, whereby 
large eddies drive all smaller ones directly, relates spa- 
tial separations r and temporal separations t linearly via 
the mean-flow velocity, whence we get trivial dynamic 
scaling with Zp = 1, for all p. A quasi-Lagrangian de- 
scription [28l [33 eliminates sweeping effects so we cal- 
culate time-dependent, quasi-Lagrangian vorticity struc- 
ture functions from our DNS. (II) Such time-dependent 
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structure functions, even for a fixed order p, do not col- 
lapse onto a scaling function, with a unique, order-p, dy- 
namic exponent. Hence, even for a fixed order p, there is 
an infinity of dynamic-multiscaling exponents [28 l [29 l [31]- 
[35j; roughly speaking, to specify the dynamics of an eddy 
of a given length scale, we require this infinity of expo- 
nents. 

Statistically steady fluid turbulence is very different in 
3D and 2D; the former exhibits a direct cascade of energy 
whereas the latter shows an inverse cascade of kinetic 
energy from the energy-injection scale to larger length 
scales and a direct cascade in which the enstrophy goes 
towards small length scales [^; in many physical real- 
izations of 2D turbulence, there is an air-drag-induced 
friction. In this direct-cascade regime, velocity structure 
functions show simple scaling but their vorticity counter- 
parts exhibit multiscaling [38-40 , with exponents that 
depend on the friction. Time-dependent structure func- 
tions have not been studied in 2D fluid turbulence; the 
elucidation of the dynamic multiscaling of these structure 
functions, which we present here, is an important step in 
the systematization of such multiscaling in turbulence. 

We numerically solve the forced, incompressible, 2D 
Navier-Stokes (2DNS) equation with air-drag-induced 
friction, in the vorticity(co')-stream-function(?/^) represen- 
tation with periodic boundary conditions: 



(1) 



where V^V^ = cj, J(V^,a;) = - (S^cj)(S^V^), 

and the velocity u = {—dyip^dx^ip)- The coefficient of fric- 
tion is II and / is the external force. We work with both 
Eulerian and quasi-Lagrangian fields. The latter are de- 
fined with respect to a Lagrangian particle, which was at 
the point at time to, and is at the position ^(^ICo, ^o) at 
time t, such that (^^(^ICo, ^o)/<^^ = u[^(t|^07 ^o), ^l, where 
u is the Eulerian velocity. The quasi-Lagrangian velocity 
field u^^ is defined 136 as follows: 



uQL(x, t|Co, to) = u[x + |(t|Co, to), t]; 



(2) 



likewise, we can define the quasi-Lagrangian vorticity 
field uj^^ in terms of the Eulerian uj. To obtain this 
quasi-Lagrangian field we use an algorithm developed in 
Ref. ^J, described briefly in the Supplementary Mate- 
rial. 

To integrate the Navier-Stokes equations we use a 
pseudo-spectral method with the 2/3 rule for the removal 
of aliasing errors [40 and a second-order Runge-Kutta 
scheme for time marching with a time step 5t = 10~^. 
We force the fluid deterministically on the second shell 
in Fourier space. And we use ja = 0.1, u = 10~^, 
and N = 2048^ collocation points [42 We obtain a 
turbulent but statistically steady state with a Taylor 
microscale A 0.2, Taylor-microscale Reynolds num- 
ber Rex 1400, and a box-size eddy-turn-over time 
'?"eddy — 8. We rcmovc the effects of transients by dis- 
carding data upto time < 80reddy We then obtain data 



for averages of time-dependent structure functions for a 
duration of time lOOreddy The energy spectrum aver- 
aged over the same time interval is shown in Fig. ([l^). 

The equal-time, order-p, vorticity structure functions 
we consider are S^{r) = {[Suj^{r^t)]P) ~ r^p , for rj^ <C 
r <^ L, where ^cj^(r,t) = [cj^(x + r,t) -cc;^(x,t)], the an- 
gular brackets denote an average over the nonequilibrium 
statistically steady state of the turbulent fluid, and the 
superscript is either E, in the Eulerian case, or QL, in 
the quasi-Lagrangian case; for notational convenience we 
do not include a subscript uj on and the multiscaling 
exponent We assume isotropy here, but show below 
how to extract the isotropic parts of 5*^ in a DNS. We 
also use the time-dependent, order-p vorticity structure 
functions 

J^^{r, {h, . . .,tp}) = {[5Lo^{r, h)... 6io^{r,tp)]); (3) 

here ti, . . . , are p different times; clearly, T^{r^ {ti = 
... = tp = 0}) = S^{r). We concentrate on the case 
ti = t2 = . . . = U = t and ti^i = ti^2 = . . . = tp = 0, 
with / < and, for simplicity, denote the resulting time- 
dependent structure function as J^^{r^t); shell- model 
studies [3T1|32] have shown that the index / does not affect 
dynamic-multiscaling exponents, so we suppress it hence- 
forth. Given J^^(r, t), it is possible to extract a charac- 
teristic time scale rp{r) in many different ways. These 
time scales can, in turn, be used to extract the order- 
p dynamic-multiscaling exponents Zp via the dynamic- 
multiscaling Ansatz rp{r) ~ r^p . If we obtain the order-p, 
degree-M, integral time scale 
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S^ir) Jo 



(l/M) 



(4) 



we can use it to extract the integral dynamic-multiscaling 
exponent z^'f^ from the relation T^'^ ^ r^vM ^ Similarly, 



from the order-p, degree-M, derivative time scale 



d 



M 



t=OJ 



(-1/M) 



(5) 



we obtain the derivative dynamic-multiscaling exponent 
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via the relation 7^ 
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Equal-time vorticity structure functions in 2D fluid 
turbulence with friction exhibit multiscaling in the di- 
rect cascade range [38ti4Q] . For the case of 3D homoge- 
neous, isotropic fluid turbulence, a generalization of the 
multifractal model [27], which includes time-dependent 
velocity structure functions [28l [32l |34l |35] , yields linear 
bridge relations between the dynamic-multiscaling expo- 
nents and their equal-time counterparts. For the direct- 
cascade regime in our study, we replace velocity struc- 
ture functions by vorticity structure functions and thus 
obtain the following bridge relations for time-dependent 
vorticity structure functions in 2D fluid turbulence with 
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friction: 

<t = l + [CtM-Cp1/M; (6) 

^^i = ^ + [(^t-CM]IM. (7) 

The vorticity field uo^ = {uo^) + uj'^ can be decom- 
posed into the time- averaged mean flow {uo^) and the 
fluctuations uo'^ about it. To obtain good statistics for 
vorticity structure functions it is important to elimi- 
nate any anisotropy in the flow by subtracting out the 
mean flow from the field. Therefore, we redefine the 
order-p, equal-time structure function to be 6'^(rc,R) = 
(|cj'^(rc + R) - cj'^(rc)l^), where R has magnitude R 
and Tc is an origin. We next use 5'p(R) = (5'^(rc, R))rc7 
where the subscript Yc denotes an average over the origin 
(we use Yc = (z,j),2 < z,j < 5). These averaged struc- 
ture functions are isotropic, to a good approximation for 
small i?, as can be seen from the illustrative pseudocolor 
plot of 5'^^(R) in Fig. ([l^). The purely isotropic parts 
of such structure functions can be obtained @0]|42 via an 
integration over the angle that R makes with the x axis, 
i.e., we calculate S^{R) = J^^ S^{lV)dO and thence the 
equal-time multiscaling exponent , the slopes of the 
scaling ranges of log- log plots of S^{R) versus R. The 
mean of the local slopes = d{\og S^) / d{\og R) in the 
scaling range yields the equal-time exponents; and their 
standard deviations give the error bars. The equal-time 
vorticity multiscaling exponents, with 1 < p < 6, are 
given for Eulerian and quasi-Lagrangian cases in columns 
2 and 3, respectively, of Table 1; they are equal, within 
error bars, as can be seen most easily from their plots 
versus p in Fig.([l]3). 

We obtain the isotropic part of J^^{R^t) in a similar 
manner. Equations Q and ([5| now yield the order- 
p, degree-M integral and derivative time scales (see the 
Supplementary Material). Slopes of linear scaling ranges 
of log- log plots of Tp'^{R) versus R yield the dynamic 

multiscaling exponent A representative plot for the 

quasi-Lagrangian case, p = 2 and M = 1, is given in Fig. 
([l]d); we fit over the range — 1.2 < \ogiQ{r/L) < —0.55 
and obtain the local slopes Xp with successive, nonover- 
lapping sets of 3 points each. The mean values of these 
slopes yield our dynamic-multiscaling exponents (column 
5 in Table 1) and their standard deviations yield the error 
bars. We calculate the degree-M, order-p derivative time 
exponents by using a sixth-order, finite-difference scheme 
to obtain T^j{f and thence the dynamic-multiscaling ex- 
ponents z^'^. Our results for the quasi-Lagrangian case 
with M = 2 are given in column 7 of Table 1. We find, 
furthermore, that both the integral and derivative bridge 
relations ([6| and ([t]) hold within our error bars, as shown 
for the representative values of p and M considered in Ta- 
ble 1 (compare columns 4 and 5 for the integral relation 



and columns 6 and 7 for the derivative relation). Note 
also that the values of the integral and the derivative 
dynamic-multiscaling exponents are markedly different 
from each other (compare columns 5 and 7 of Table 1). 

The Eulerian structure functions Tp{R,t) also lead 
to nontrivial dynamic-multiscaling exponents, which are 
equal to their quasi-Lagrangian counterparts (see Sup- 
plementary Material). The reason for this initially sur- 
prising result is that, in 2D turbulence, the friction con- 
trols the size of the largest vortices, provides an infra-red 
cut-off at large length scales, and thus suppresses the 
sweeping effect. We have demonstrated this in the sup- 
plementary material. Had the sweeping effect not been 
suppressed, we would have obtained trivial dynamic scal- 
ing for the Eulerian case. 

The calculation of dynamic-multiscaling exponents has 
been limited so far to shell models for 3D, homogeneous, 
isotropic fluid |29l [31] [32j |^ 35 and passive-scalar tur- 
bulence [33]. We have presented the first study of such 
dynamic multiscaling in the direct-cascade regime of 2D 
fluid turbulence with friction by calculating both quasi- 
Lagrangian and Eulerian structure functions. Our work 
brings out clearly the need for an infinity of time scales 
and associated exponents to characterize such multiscal- 
ing; and it verifies, within the accuracy of our numerical 
calculations, the linear bridge relations (|6| and ^ for a 
representative value of /i. We find that friction also sup- 
presses sweeping effects so, with such friction, even Eu- 
lerian vorticity structure functions exhibit dynamic mul- 
tiscaling with exponents that are consistent with their 
quasi-Lagrangian counterparts. 

Experimental studies of Lagrangian quantities in tur- 
bulence have been increasing steadily over the past 
decade [44 . We hope, therefore, that our work will en- 
courage studies of dynamic multiscaling in turbulence. 
Furthermore, it will be interesting to check whether the 
time scales considered here can be related to the persis- 
tence time scales for 2D turbulence [45 . 
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FIG. 1. (Color online) (a) Log-log (base 10) plot of the energy spectrum E{k) versus k. The black line with slope —4.1 is 
shown for reference, (b) Pseudocolor plot of the equal-time, quasi-Lagrangian, second-order vorticity structure function *S'^^(R) 
averaged over the origin rc (we use rc = {i,j),2 < i,j < 5); (c) plots of the equal-time, vorticity, multiscaling exponents 
versus p for Eulerian (red circles) and quasi-Lagrangian (blue diamonds) fields (error bars are comparable to the size of the 
symbols); the inset shows the local slopes ^p, obtained as defined in the text, versus the separation, from p = 1 (bottom) to p 
— 6 (top); (d) log-log (base 10) plot of the order-2, degree-1, integral time scale T2 '^^{R) versus the separation R showing our 
data points (open red circles) and the best-fit line (full black) in the scaling range; the inset shows the local slopes Xp^ obtained 
as defined in the text, versus the separation, from p = 1 (bottom) to p = 6 (top). 



order (p) 




Sp 




pA 




D,QL 


1 

2 
3 
4 
5 
6 


0.62 zb 0.009 
1.13 it 0.009 
1.561 ± 0.009 
1.92 ± 0.01 
2.24 ± 0.01 
2.52 ± 0.02 


0.63 ± 0.008 
1.13 ± 0.008 
1.54 ± 0.01 
1.90 ± 0.01 
2.26 ± 0.01 
2.54 ± 0.02 


0.366 ± 0.008 
0.50 ± 0.02 
0.59 ± 0.02 
0.64 ± 0.02 
0.64 ± 0.02 
0.72 ± 0.03 


0.37 =b 0.02 
0.48 ± 0.01 
0.57 ± 0.01 
0.63 ± 0.01 
0.65 ± 0.02 
0.67 =b 0.02 


0.55 =b 0.02 
0.62 ± 0.02 
0.64 ± 0.02 
0.68 =b 0.03 
0.70 =b 0.03 
0.71 ± 0.04 


0.53 ± 0.02 
0.62 ± 0.2 
0.65 ± 0.01 
0.68 ± 0.01 
0.70 ± 0.02 
0.71 ± 0.03 



TABLE L Order-p (column 1); equal-time, Eulerian exponents Cp (column 2); equal-time, quasi-Lagrangian exponents 
(column 3); integral-scale, dynamic- multiscaling exponent z^^^^ (column 4) from the bridge relation and the values of in 
column 3; z^^^^ from our calculations of time-dependent structure functions (column 5); the derivative-time exponents z^J^^ 
(column 6) from the bridge relation and the values of in column 3; z^J^^ from our calculations of the time-dependent 
structure function (column 7). The error estimates are obtained as described in the text. 
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Algorithm for obtaining quasi-Lagrangian fields in a 
pseudospectral simulation 

To obtain a quasi-Lagrangian field from its Eulerian 
counterpart, we track a single Lagrangian particle by us- 
ing a bilinear-interpolation method [17]. If we replace 
the Eulerian velocity in Eq. (2) by its Fourier-integral 
representation, we obtain 

uQL(x,t|Co,^o) = j u(q,t)exp[iq-(x + C(t|Co,^o),^))]^^q, 

where q is the wave vector. In the pseudospectral al- 
gorithm we use to solve Eq. (1), the quasi-Lagrangian 
velocity is defined with respect to a Lagrangian parti- 
cle, which was at the point at time to, and is at the 
position C(t|^07^o) at time t, such that d$^{t\$^i^^tQ) / dt = 
u[^(t|^07 ^o), We calculate u(q, t); thus, the Fourier 
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integral above can be evaluated at each time step by an 
additional call to a fast-Fourier- Transform (FFT) sub- 
routine. The additional computational cost of obtaining 
u^^ at all collocation points is that of following a single 
Lagrangian particle and an additional FFT at each time 
step. 



Numerical determination of integral time scales 
from time-dependent structure functions 

To extract the integral time scale, of degree M, from a 
time-dependent structure function, we have to evaluate 
the integral in Eq. (4) numerically. In practice, because 
of poor statistics at long times, we integrate from t = 
to t = t*, where is the time at which T^{R\t) = e; 
we choose e = 0.6, but we have checked that our results 
do not change, within our error bars, for 0.5 < e < 0.75. 
This numerical integration is done by using the trape- 
zoidal rule. 



Dynamic multiscaling for Eulerian structure 
functions 

Equal-time Eulerian structure functions have been dis- 
cussed in our paper above. To obtain time-dependent, 
Eulerian, vorticity structure functions we proceed as we 
did in the quasi-Lagrangian case. We obtain the required 
vorticity increments and from these the purely isotropic 
part of the time-dependent, order-p structure function 
Tp{R'^t). Equations (4) and (5) now yield the order- 
p, degree-M integral and derivative Eulerian time scales. 
For the former we should integrate Tp{R'^t) from t = 
to t = oo; in practice, because of poor statistics at long 
times, we integrate from t = to t = t*, where is 
the time at which J^^{R' ^t) = e; we choose e = 0.6, but 
we have checked that our results do not change, within 
our error bars, for 0.5 < e < 0.75. Slopes of linear scaling 
ranges of log- log plots of T^'^{R') versus R' yield the dy- 
namic multiscaling exponent z^'f . A representative plot 
for the Eulerian case, p = 2, and M = 1 is given in Fig. 
([2] a); we fit over the range —1.2 < log;Lo(^/^) < —0.55 



and obtain the local slopes Xp with successive, non- 
overlapping sets of 3 points each. The mean values of 
these slopes yield our dynamic-multiscaling exponents 
(column 4 in Table [ll|) and their standard deviations 
yield the error bars. We calculate the degree-M, order- 
p derivative time exponents by using a sixth-order fi- 
nite difference scheme to obtain and thence the 

D E 

dynamic-multiscaling exponents z^Jj^] data for the Eule- 
rian case and the representative value M = 2 are given 
in column 6 of Table [HI We find, furthermore, that both 
the integral and derivative bridge relations, Eq. (6), and 
Eq. (7). hold within our error bars, as shown for the 
representative values of p and M considered in Table |n| 
(compare columns 3 and 4 for the integral relation and 
columns 5 and 6 for the derivative relation). The values 
of the integral and the derivative dynamic-multiscaling 
exponents are markedly different from each other (com- 
pare columns 4 and 6 of Table |ll| and the plots of these 
exponents versus p in Fig. ([2]b). In Fig. ([2]c), we make 
the same comparison for the quasi-Lagrangian case. Fur- 
thermore, a comparison of the quasi-Lagrangian and Eu- 
lerian dynamic-multiscaling exponents given in Tables I 
in the original paper and Table [Til respectively, show that 
these are the same (within our error bars). 

Demonstration of Infra-red cutoff of the inverse 
cascade 

We have shown that in two dimensional turbulence 
with friction, the Eulerian and the quasi-Lagrangian ve- 
locities have the same dynamical exponents. This is be- 
cause the inverse cascade has a friction-dependent infra- 
red cutoff. 

To illustrate the development of this cutoff scale, we 
have carried out DNS studies of 2D fluid turbulence with 
/i = 0.01, 0.05, and 0.1, 1024^ collocation points, and 
forcing at a wave- vector magnitude /c = 80; our DNS 
studies resolve the inverse-cascade regime in the statisti- 
cally steady state. The energy spectra from these DNS 
studies, plotted in Fig. ([2]i), show clearly that, as /i in- 
creases, the inverse cascade is cut off at ever larger values 
of k. Thus, the friction produces a regularization of the 
flow and suppresses infrared (sweeping) divergences. 
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FIG. 2. (Color online) (a) Log-log (base 10) plot of the order-2, degree-1, integral time scale tI'^[R') versus the separation 
R' showing our data points (open red circles) and the best-fit line (full black) in the scaling range; the inset shows the local 
slopes obtained as defined in the text, versus the separation, from p — 1 (bottom) tp p = 6 (top); (b) plots of the vorticity, 
dynamic- multiscaling, quasi- Lagrangian exponents 2;^'^^ (open red circles) and z^ J^^ (full blue circles) versus p with the error 
bars given in columns 4 and 6, respectively, in Table I; (c) plots of the vorticity, dynamic- multiscaling, Eulerian exponents z^'f 
(open red circles) and z^J^ (full blue circles) versus p with the error bars given in columns 4 and 6, respectively, in Table 11. 
(d) log-log (base 10) plot of the energy spectrum E(k) versus the wave- vector magnitude k for /x = 0.01 (red filled circles), 
/i = 0.05 (blue filled diamonds), and /i = 0.1 (magenta filled triangles); the peak is at the injection scale /c = 80 and the black 
line indicates the K41, 2D-inverse-cascade slope. 



order (p) 




<f[Eq.(6)] 


/,E 


<f[Eq-(7)] 




1 


0.62 ± 0.009 


0.372 ± 0.009 


0.37 =b 0.02 


0.53 zb 0.02 


0.53 ± 0.02 


2 


1.13 ± 0.009 


0.51 ± 0.02 


0.48 ± 0.01 


0.60 ± 0.02 


0.62 zb 0.02 


3 


1.561 ± 0.009 


0.56 ± 0.02 


0.57 zb 0.01 


0.66 lb 0.02 


0.65 zb 0.02 


4 


1.92 ± 0.01 


0.64 ± 0.02 


0.63 zb 0.01 


0.70 zb 0.03 


0.68 zb 0.02 


5 


2.24 ± 0.01 


0.68 ± 0.02 


0.65 ± 0.02 


0.71 zb 0.03 


0.70 zb 0.02 


6 


2.52 ± 0.02 


0.72 ± 0.03 


0.67 ± 0.02 


0.72 ± 0.03 


0.71 ± 0.03 



TABLE IL Order-p (column 1); equal-time, Eulerian exponents (column 2); integral-scale, dynamic- multiscaling exponent 
Zp'f (column 3) from the bridge relation and the values of ("p in column 2; ^^'f from our calculation of time-dependent structure 
functions (column 4); the derivative-time exponents z^J^ (column 5) from the bridge relation and the values of (p in column 
2; z^2^ from our calculation of time-dependent structure function (column 6). The error estimates are obtained as described 
in the text. 



